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ABSTRACT 


The  supersonic  flow  field  about  a  forebody- 
decelerator  combination  is  analyzed.  Tangent-cone  and 
Newtonian  methods  are  used  to  predict  forebody  surface 
pressures.  Empirical  correlations  based  on  local 
similarity  solutions  predict  the  boundary  layer  character¬ 
istics,  and  boundary  layer  similarity  techniques  are  used 
to  predict  the  forebody  wake  parameters.  The  flow  about 
a  conical  secondary  body  in  the  forebody  wake  is  calculated 
using  the  method  of  characteristics.  In  Volume  I  the 
methods  are  presented  and  explained,  and  the  results  of  the 
calculations  are  compared  with  experimental  measurements; 
in  Volume  II  the  details  of  the  calculation  procedure  and 
computer  program  are  presented. 
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I.  INTRODUCTION 


This  report  presents  the  results  of  an  effort  to 
calculate  the  flow  field  characteristics  of  a  two-body  trailing 
decelerator  system  in  supersonic  flow  (Fig  1) .  The  report 
is  divided  into  two  volumes,  the  first  dealing  with  theoreti¬ 
cal  methods  and  a  comparison  of  the  results  of  a  sample  cal¬ 
culation  with  measurements,  the  second  presents  and  explains 
the  computer  program  used  for  the  calculations. 

This  volume,  Volume  I,  presents  an  explanation  of 
the  methods  used  to  analyse  the  forebody  (the  object  to  be 
decelerated  or  payload) ,  the  trailing  secondary  body  (the 
decelerator)  (Fig  1) ,  and  the  comparison  of  results  with 
experimentally  measured  values.  Naturally,  the  ideal  objective 
would  be  to  have  a  theory  and  methods  that  provide  a  complete 
closed  solution  to  the  entire  flow  field  about  the  two  bodies. 
Within  the  curre'-*  state  of  the  art,  however,  this  goal  cannot 
be  achieved.  Th.  >  . ^forts  have  been  directed  at  utilizing 
tractable  method  *dth  limited  assumptions  and  required  inputs 
that  provide  very  useful  results,  wherever  assumptions  are 
necessary,  they  are  discussed  in  detail  and  their  effects 
evaluated  where  possible. 
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II.  FOREBODY  AND  WAKE  ANALYSIS 


This  division  describes  the  analysis  used  to 
predict  the  flow  characteristics  about  the  forebody  and  its 
wake,  and  draws  heavily  on  the  work  presented  in  Ref  1,  in 
particular  Section  B,  the  explanation  of  the  analysis. 

In  view  of  this,  the  matter  presented  in  this  section  is  not, 
nor  is  it  intended  to  be,  rigorously  complete;  rather  it 
should  prov5de  a  basic  understanding  of  the  analysis.  For 
further  details  or  information  Ref  1  must  be  consulted. 

A.  Extent  of  Analysis,  Required  Inputs,  and  Results 

1.  Extent  of  Analysis 

The  ^orebody  and  wake  analysis  examines  the  flow 
about  the  forebody  extending  axially  from  the  bow  shock  down¬ 
stream  to  the  secondary  body  bow  shock,  and  normally  from  the 
surface  of  the  forebody  outward  to  its  bow  shock.  The  analysis 
assumes  tha^  flow  conditions  upstream  of  the  forebody  are 
uniform. 

2.  Required  Inputs 

The  forebody  and  wake  analysis  requires  the  following: 

a.  freestream  conditions 

b.  forebody  geometry 

c.  bow  shock  shape  and  standoff  distance. 

3.  Results 

In  addition  to  body  geometry  and  freestream  condi¬ 
tions,  the  forabody  subroutine  produ,*.s  the  calculated  local 
static  pressure,  total  pressure  and  Mach  number  on  the  body 
surface  at  various  axial  locations  ar  i  on  various  planes  con¬ 
taining,  and  rotated  about,  the  sxis.  Since  the  flow  char¬ 
acteristics  at  the  base  of  the  fotebody  are  likely  to  be 
affected  by  a  possible  transition  into  turbulence  on  the  body 
surface,  the  transition  point  location  on  the  body  surface 
and  on  various  planes  is  also  printed  out.  Reynold's 
number,  based  on  momentum  thickness,  inviscid  density,  local 
inviscid  velocity,  local  momentum  defect,  and  momentum  thick¬ 
ness  at  the  base  for  different  $  planes,  is  provided.  Inte¬ 
gration  of  the  local  momentum  defect  around  the  surface  at 
the  forebody  base  provides  the  total  momentum  defect,  which 
will  be  used  in  wake  computations. 


Wake  computations  begin  with  forebody  boundary  layer 
and  base  condition  inputs.  The  program  selects  either 
laminar  or  turbulent  wake  conditions,  based  on  the  transition 
Reynold's  number,  and  prints  out  the  input  values  of  the 
viscosity  model.  At  every  JU D  location  behind  the  base  of 
the  forebody  and  at  every  iteration  of  wake  computation, 
values  of  ©/A,  ,  SI D,  <f  *,  and  the  local  inviscid  velocity 

at  the  wake  edge  are  provided.  The  viscous  wake  is  specified 
by  the  values  of  R/D,  M,  M/M„  ,  V/V  ,  u,  p  /p^,  h  Jh  ,  and 

-  U  O  Oc 

h/h^  at  each  axial  location,  ^/D,  and  radial  distance  ^  ,  the 

non-dimensional  transformed  radial  coordinate.  The  inviscid 
wake  is  determined  by  the  values  of  R/D,  M,  M/M® 
u,  P0/p®  >  P0/Poe*  anc*  T  at  var*ous  values  of  the 

coordinate  yt.  . 


,  v/v 

transformed 


B.  Method  and  Explanation  of  the  Forebody-Wake 

Analysis  (from  Ref  1) 

1.  Body  to  Bow  Shock  Region 


The  variations  of  flow  properties  in  this  region  are 
primarily  due  to  the  changes  caused  by  the  forebody  bow  shock 
and  secondary  shocks  such  as  recompression  shocks  and  forebody 
flare  shocks.  Assuming  negligible  changes  in  entropy  across 
the  secondary  shocks,  total  flow  properties  are  invariant 
following  a  streamline  starting  from  any  point  behind  the 
forebody  bow  shock.  In  Ref  1,  the  shock  shape,  determined  from 
experiments  or  from  analytical  methods  for  simple  bodies,  is 
represented  in  multiple  regression  curve  fit  form 

r_  Q  *  b(lk)  *  °(%'s)  (1) 

'-d(D 

where  a,  b,  c,  and  d  are  polynomial  coefficients.  Since  the 
shock  is  completely  specified,  the  flow  properties  immediately 
behind  the  shock  are  determined  from  shock  jump  conditions. 

2.  Forebody  Surface  Pressure  Distribution 


The  total  and  static  pressures  on  the  forebody  sur- 
face  are  estimated  using  Newtonian  or  tangent-cone  methods. 

The  forebody  surface  is  represented  by  a  number  of  x  and  r 
coordinates  of  points  on  the  surface.  The  nose  half- angle  of  the 
forebody  is  determined  using  coordinates  at  Oth  and  1st 
stations  where 

<rL  —  tan'1  (—*  _ 1— j  at  t  *0  (2) 
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For  a  given  freestream  Mach  number  an  attached  shock  exists 
if  the  nose  half-angle  is  less  than  the  critical  cone  half¬ 
angle.  The  relation  between  critical  cone  half-angle  and 
freestream  Mach  number  (Ref  2)  can  be  represented,  after  curve 
fitting,  by 

^  I  -  *  HMH  rt3  -  /.  07148  M4 

0C  =  - - -  (3) 

(I-  1 11017 M) 

If  the  nose  half-angle  £  from  Eqn  2  is  greater  than  or 
equal  to  the  critical  cone  angle  from  Eqn  3,  a  detached  shock 
exists  at  the  apex.  The  apex  total  pressure  p  is  given  by 
normal  shock  relation  (Ref  2) 


& 


J(MA l]*  r  y+i  i  ) 

(*•  1  j  IztMf  -Of-oJ  ( 


(4) 


Static  pressure  can  be  estimated  from  the  Newtonian  method, 


&  ~  yf  Sin2^  (5) 

If  the  nose  half-angle  <T0  is  less  than  critical  cone  half¬ 
angle  £  ,  an  attached  shock  exists  at  the  apex.  The  flow 
region  behind  the  shock  is  considered  a  conical  flow  region. 
The  nose  half-angle  $  and  the  freestream  Mach  number  approx¬ 
imately  determine  the  shock  angles  .  The  velocity  components 
along  and  perpendicular  to  any  ray  from  the  apex  of  the  body 
can  be  determined  by  expanding  the  velocities  in  Taylor  series 
and  using  the  approximate  shock  angle  ft.  The  procedure  is 
continued  until  the  normal  velocity  on  the  body  surface  is 
evaluated.  Since  the  boundary  condition  at  the  surface 
requires  the  normal  component  of  velocity  to  be  zero,  the 
shock  angle  /3  is  altered  until  the  deviation  of  normal 
velocity  from  zero  is  within  allowable  limits.  The  surface 
Mach  number,  static  pressure,  and  total  pressure  are  evaluated 
from  the  above  conical  flow  analysis. 

A  Newtonian  or  conical  flow  solution  is  used  to 
determine  the  static  pressure  at  each  point  on  the  body  pro¬ 
file.  As  long  as  the  local  inclination  angle  dV*,  of  the 
body  surface  at  the  (i  +  l)th  station  is  less  than  or  equal  to 
the  local  inclination  angle  tf}  at  the  ifch  station,  no  surface 
shock  exists  and  total  pressure  remains  unaltered.  If  the 
bow  shock  wave  is  detached,  static  pressure  is  obtained  from 
the  Newtonian  method,  until  the  local  inclination  angle  62  is 
less  than  the  critical  cone  half-angle  6^  beyond  which  the 
conical  flow  solution  is  used. 


On  the  other  hand,  if  the  local  inclination  angle 
(Ti+,  at  the  (i  +  l)t“  station  is  greater  than  the  value  <57 
at  the  ith  station,  a  corner  shock  appears  at  the  surface, 
and  the  total  pressure  is  altered.  Tne  effective  conical 
half- angle,  5=  ,  is  given  by 


=  <i.,  -  c 

Here  again,  if  W  is  greater  than  or  equal  to  the  critical 
half-angle  corresponding  to  Mach  number  Mgi  at  the  i™ 

station,  the  shock  is  detached,  and  the  total  pressure 


(6) 


cone 


P0(i  +  1)  at  (i  4  1)^  station  will  be 


=  ft, 


(i*0 


1* 

[  M  l] 

[(H)  4,/ ->2. 

M 


(7) 


The  static  pressure  is  determined  from  the  Newtonian  method, 


ft0  Sml  <T 
(Un 


(8) 


However,  if  ?  is  less  than  critical  cone  half-angle  <rc 
corresponding  to  Me^,  the  conical  flow  solution  will  be  used 

to  determine  total  pressure,  static  pressure,  and  surface 
Mach  number. 


3.  Boundary  Layer  Analysis 

Because  of  the  nature  of  the  wake  analysis,  IIB4, 
the  primary  aim  of  the  boundary  layer  analysis  is  to  obtain 
the  momentum  defect  at  the  base  of  the  forebody,  which  is 
also  the  only  output  from  the  computer  program.  Thus  this 
boundary  layer  analysis  is  restricted  to  determining  that 
quantity,  and  assumes  that  the  flow  properties  at  the  edge 
of  the  boundary  layer  are  those  derived  in  1IB2  for  forebody 
surface  conditions. 

The  boundary  layer  momentum  defect  can  be  expressed 
as 


m  -  %  ue 


/■>  /* 


O  • O 


(9) 


Equation  9  can  be  recast  into  the  form 

m*  -  j'  fe  O  ds  (io) 

Jarc  u  m ferencz 

where  ds  is  the  circumferential  increment  of  length.  The 
boundary  layer  momentum  thickness  e  is  defined  as 


Since  boundary  layer  characteristics  primarily  depend  on  the 
laminar  or  turbulent  nature  of  the  flow  on  the  forebody  sur¬ 
face,  the  momentum  thickness,  e  ,  for  these  two  cases  is 
developed  in  the  following  paragraphs,  a  and  b,  and  the  laminar- 
turbulent  transition  is  discussed  in  paragraph  c. 

a.  Laminar  Boundary  Layer  Analysis 

The  momentum  thickness  of  a  laminar  forebody  boundary 
layer  is  determined  from  an  empirical  equation  which  depends 
on  edge  conditions  and  local  pressure  gradient. 

This  effect  of  pressure  gradient  on  a  laminar 
boundary  layer  is  usually  taken  in  the  form  of  the  pressure 
gradient  parameter/?  ,  which  varies  from  one,  at  a  two- 
dimensional  stagnation  point,  to  zero  for  a  zero  pressure 
gradient  region.  Although  a  varies  from  point  to  point  on 
the  body  surface,  the  bounaary  layer  characteristics  largely 
depend  on  the  local  value  of  the  pressure  gradient  parameter fi  . 
This  concept  of  "local  similarity"  has  been  used  in  Ref  1  to 
obtain  the  following  expression  for  e  : 


©  vr(am(t-  0.0  9/s  a4)  5*  1/2 _ 

'S  ~  Ifc  Me  /  M<  \'4  )  'kir  \K  I  ft /“C  f"4  (12> 

I  /“  /  V  /  \-Ps  I  [rel  ( I 

where  K  «  0  for  two-dimensional  flow,  and  K  1  for  axi- 
symmetric  flow. 


The  pressure  gradient  parameter  ft  is  given  by 


(13) 
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J^r 

represents  an  effective  surface  distance  considering 
boundary  layer  history  up  to  that  point, 


Knowing  the  inviscid  flow  properties  at  the  edge  of  the 
boundary  layer  and  evaluating  J*t  /?  and  0/t6  from  Eqns  14, 
13,  and  12,  the  momentum  defect  at  the  base  of  the  circular 
body  can  be  written  as 

* '  *  *T  Da{  fs  Uc  &)B  (15) 


b.  Turbulent  Boundary  Layer  Analysis 

The  turbulent  boundary  layer  analysis  used  in  Ref  1 
to  determine  the  momentum  defect  at  the  base  of  the  forebody 
is  mainly  based  on  the  results  of  Reshotko  and  Tucker  (Ref  3) . 
If  the  flow  on  the  forebody  surface  is  completely  turbulent, 
the  turbulent  boundary  layer  analysis  is  applied  from  the  nose 
stagnation  point.  However  if  the  boundary  layer  is  initially 
laminar,  turbulent  flow  analysis  starts  from  the  transition 
point  where  0  turbulent  «  ©  laminar.  The  transition  point 
is  located  on  the  basis  of  transition  Reynold's  number, 

Re®TR*  discussed  in  Section  c.  In  Ref  3  the  presence  of 

eddy  viscosity  in  the  turbulent  boundary  layer  necessitates 
the  use  of  an  approximate  equation  to  represent  the  shear 
stress  distribution,  and  the  effect  of  pressure  gradient  on 
boundary  layer  velocity  profile  is  taken  into  account  through 
the  use  of  moment  of  momentum  equation.  After  considerable 
simplification,  the  results  of  Reshotko  and  Tucker  (Ref  3 ) 
are  represented  by  a  transformed  momentum  thickness  &  tr  in 

the  mathematical  form  (Ref  1) 
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The  subscript  TR  indicates  the  values  at  the  transition  point; 
a0  and  ^  are  the  velocity  of  sound  and  the  kinematic  vis¬ 
cosity  at  stagnation  conditions  of  the  local  external  stream. 
The  turbulent  distance  parameter  j*  is  given  by 


The  distance  s^  is  the  location  of  the  transition  point 
measured  along  the  body  surface  from  the  front  stagnation 
point.  If  the  flow  is  completely  turbulent  si  *  0.  The 
exponent  B  in  Eqns  16  and  17  is  approximated  by  the  function: 


B  =  4.i  *  8  [  b  - 1  j  o 

where  TQ  is  the  inviscid  stagnation  temperature.  The  term 
ijnf  J  in  Eqn  17  is  given  as 

(fy-  " 

The  actual  momentum  thickness  Q  is  then  evaluated  from  the 
relation 

*m  9»0  *  <2 


The  local  momentum  thickness  9  for  a  turbulent  boundary  layer 
beyond  the  transition  point  can  then  be  determined  using  the 
above  equations.  The  momentum  defect  at  the  base  of  the  fore¬ 
body  is  then  determined  from  Eqn  15. 

c.  Boundary  Layer  Transition 

If  the  flow  on  body  surface  is  initially  laminar, 
the  transition  point  must  be  located  to  change  the  calculation 
method  from  the  laminar  to  the  turbulent  boundary  layer 
analysis.  The  Reynold's  number  based  on  momentum  thickness 
at  the  transition  point  can  be  approximated  by  (Ref  1) 


9 


Ree  Z  toootl-aWM1  +  Q.0ZSH3)ll%['Zl9  + 

**  '  •  \2  ,i  mi  h. 


+  17-38 


ft) 


-I8i 


k 


+  3-91 


9oi  , 


ft) 


(21) 


where  hw  =  enthalpy  based  on  wall  temperature  and 
enthalpy  based  on  adiabatic  wall  temperature. 


aw 


In  an  adverse  pressure  gradient  region,  Re0TR  Is 

replaced  by  the  critical  Reynold’s  number,  Re/*  .  given 
in  Ref  1  cr 


a  lb')  R& 


(22) 


Equations  21  and  22  have  been  developed  and  discussed  in 
greater  detail  in  Appendix  III  of  Ref  1.  The  transition 
point  is  located  where  the  laminar  $  satisfies  Eqns  21  or  22. 


4.  Wake  Analysis 


a.  Viscous  Wake 


The  viscous  wake  can  be  divided  into  three  major 
regions  1)  the  near  wake,  2)  the  neck  region  and  rear  stag¬ 
nation  point,  and  3)  the  far  wake.  The  near  wake  is  char¬ 
acterized  by  l)  outer  rotational  inviscid  flow,  2)  a  free 
shear  layer  above  the  dividing  stream  line  and  3)  a  recircula¬ 
tion  region.  In  the  near  wake  and  wake  neck  region,  because 
of  the  presence  of  a  pressure  gradient,  the  existence  of  a 
rear  shock  and  the  nonsimilar  nature  of  the  governing  equations, 
not  much  success  could  be  achieved  in  representing  the  velo¬ 
city  and  enthalpy  profiles  in  a  mathematical  form.  In  the 
far  wake  region  assumptions  such  as  constant  momentum  defect, 
negligible  pressure  gradient,  and  boundary  layer  assumptions 
could  be  used  to  determine  the  flow  variables.  Neglecting  non- 
isentropic  processes  in  the  near  wake  region  and  assuming 
that  no  external  work  is  done  on  the  fluid,  and  a  constant  wake 
momentum  defect,  the  wake  momentum  defect  may  be  equated  to 
the  boundary  layer  momentum  defect  at  the  base  of  the  forebody. 
Thus,  momentum  thickness,  ,  of  the  wake  is  given  by 
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(23) 


In  order  to  study  the  flow  field,  the  governing  differential 
equations  for  a  laminar  axisymmetric  high-speed  wake  are 
transformed  from  physical  coordinates  (x,r)  to  Dorodnitzyn 
coordinates  (x,w  ),  which  are  related  by 


■7 

Ldl 


'Tj  J  V  - 


±_ 

A1 


where 


A*  = 


I—  r  Jr 


J—  r  6  r 

f. 


(24) 


(25) 


The  velocity  prof ile  in  the  viscous  wake  can  be  assumed  to  be 
a  fourth  order  polynomial  in^  ,  as 

M.  =  Q0  *  Q,^  -*  <2j^5+  a,  (26) 


To  determine  the  coefficients  in  Eqn  26,  the  following 
boundary  conditions  may  be  imposed: 


At  the  centerline  ft  » 


0 


(1)  (A  - 

(2)  =  o  (27) 

(3)  axial  direction  momentum 
equation  applies 


At  the  outer  edge  of  the 
wake  ^  «  1 


(4) 

(5) 

(6) 


U*  U, 

-  o 

yv 


*  o 


(28) 


The  above  six  boundary  conditions  are  needed  to  determine 
six  unknowns,  namely  aQ  to  a^,  and  n  ,  provided  the  flow 

properties  at  the  outer  edge  of  viscous  wake  are  known  from 
the  inviscid  wake  solution.  Applying  the  above  boundary 
conditions,  the  velocity  profile  in  the  transformed  plane 
becomes 


U.  ~  §4  ('~  ) 


(29) 
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where 


5  =  ^  A2  1l 

Mi  /* ct  d%  (30) 

and  also 

5  =  24  ("  (31) 

IC  should  be  noted  that: 


(!)  * 

at  the  wake  stagnation  point 

X  -  0, J  =  24,  ucX  -  0 

(2) 

at  the  wake  neck  ^  -  0,  *•  22,9 

^  °  0.045 

4 At 

(32) 

(3) 

far  downstream,  X—  <*>,  $  »*  0, 

iii i  . 

uT  u  1 

In  order  to  specify  the  velocity  profile  completely  at  any 
location  from  the  base  of  forebody ,  it  is  necessary  to  know 
the  distance  of  the  wake  stagnation  point  *•  ,  from  the  nose 
of  the  forebody,  and  the  relation  between  ?  ,  the  form  para¬ 
meter,  and  X,  the  longitudinal  distance  from  the  wake  stagna¬ 
tion  point.  It  is  difficult  to  predict  the  location  of  a 
wake  stagnation  point  from  mathematical  analysis,  however, 
from  the  limited  number  of  experiments  for  turbulent  flow, 
it  is  taken,  in  Ref  1,  that  the  wake  stagnation  point  is 
located  approximately  at  0,94  Dg  from  the  base  of  the  forebody. 

The  relation  between  jT  and  X  will  then  be  obtained  from  Rqn  30. 
Because  of  the  presence  of  the  viscosity  coefficient/^ in 
Eqn  30  the  relation  between  $  and  X  depends  on  the  laminar 
or  turbulent  nature  of  the  flow  in  the  wake.  These  are  dis¬ 
cussed  in  IIB4a(l)  and  HB4a(2). 


The  assumption  of  zero  pressure  gradient  in  the 
wake  and  a  Prandtl  number  of  unity,  leads  to  a  solution  of 
the  energy  equation  in  the  form 


Uct 

U, 


(33) 


The  total  enthalpy,  at  the  wake  centerline  can  be 

determined  using  the  condition  that  the  total  enthalpy  defect 
in  the  wake  is  constant  and  is  equal  to  the  value  at  the  base 
of  the  forebody.  Carrying  out  the  algebra,  it  can  be  shown 
that 


where  tu  is  the  enthalpy  based  on  the  skin  temperature  of 
the  forebody  at  the  base.  Substituting  for  hoc^  from  Eqn  34 

and  into  Eqn  33,  we  get  the  total  enthalpy  profile 

Knowing  the  velocity  profile  from  Eqn  29  and  the  total  enthalpy 
profile  from  Eqn  35,  the  static  enthalpy  profile  can  be  ob¬ 
tained  from  the  relation 


h-  4  -  £  <36> 

z 

The  assumption  of  wake  static  pressure  equal  to 
freestream  static  pressure  and  the  velocity  and  static  enthalpy 
profiles  as  determined  from  29  and  36  are  sufficient  to  deter¬ 
mine  all  the  flow  properties  in  the  viscous  wake.  Since  the 
profiles  are  determined  in  terms  of  Dorodnitzyn  coordinate  7  , 
it  is  then  necessary  to  transform  back  into  physical  coordinates, 
using  the  relation 


=  l 


L  h* 


(1.)  Laminar  Viscous  Vake 

In  order  to  specify  the  velocity  and  total  enthalpy 

?rof iles  at  any  specified  distance  from  the  base  of  the  forebody, 
t  Is  necessary  to  determine  the  relation  between  $  *  the  form 
parameter,  and  X,  the  longitudinal  distance  from  wake  stag¬ 
nation  point.  For  a  laminar  flow  case,  a  linear  relationship 
between  viscosity  and  flow  field  temperature  can  be  assumed 
(i.e. ,  y^/RT  *  constant).  Because  of  the  validity  of  the 
boundary  layer  assumptions  in  the  viscous  wake,  we  get  from 
the  equation  of  state, 

f* 


13 


(38) 


The  density  at  the  centerline  of  wake  is  given  by 

\C  sjp 

UK 

From  the  definition  of  the  momentum  equation  in  transformed 
coordinates, 


JN  (40) 

Substituting  the  velocity  profile  from  Eqn  29  into  40,  we 
get  the  relationship  between  ^  ,  the  wake  thickness,  and  T  , 
the  form  parameter. 

a.1  ,,  y_  ,  xl  (4i) 

£*  /ZO  7?O0 

Using  Eqns  38,  39,  41,  36,  31,  and  30,  we  get  a  first  order  non- 
linear  differential  equation 

£S__JL  1 

Y-!  /*•  &  24  ‘  241  \  120  5700) 

(42) 


Integrating  Ean  42  from  X  *  0, J  -  24,  we  get  the  relationship 
between  $  ana  X  for  a  laminar  viscous  wake: 


±_  i>,e'(  ill  & Ldt  _ 

h  *♦  j/r.  1  ui-hi\-u't,-i- 

sst»\  1  us  14  hi  I  TV  3 


im 


(2.)  Turbulent  Viscous  Wake 

The  presence  of  eddy  viscosity  in  turbulent  flow 
necessitates  the  use  of  an  effective  viscosity  coefficient 
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in  Eqn  30  to  determine  the  relationship  between  $  ,  the  form 
parameter,  and  X,  the  longitudinal  distance  from  the  wake 
stagnation  point.  Four  viscosity  models  have  been  suggested 
in  Refs  1  and  4,  however,  the  third  model,  namely 

'  k  Af ,  (  Ut~Ua)  (44) 


with  a  value  of  k  «  0.005  is  found  to  give  the  best  correla¬ 
tion  with  their  experimental  results.  Substituting  Eqn  44 
into  30,  replacing  u  in  terms  of  5  >  f*om  Eqns  31  and 

41,  we  get  a  first  order  ordinary  differential  equation  in 
the  form  (for  the  third  viscosity  model) 


d_x  =  £  fimr 

dz.  ~0,U-  Juj  V /zo  55i,° 


(45) 


Assuming  the  wake  momentum  thickness,  9,  ,  to  be  a  constant, 
and  integrating  Eqn  45  from  XTR,  the  turbulent  wake  transition 

location,  we  get  the  relationship  between  $  and  X  as 


where  XTR  »  location  of  wake  transition  point  and  TR  -  £ 
laminar  at  transition  point. 


Similarly,  the  relationship  between  Jr  and  X  can  be 
obtained  for  other  viscosity  models,  but  in  a  more  complicated 
form.  If  the  flow  on  the  forebody  surface  is  turbulent,  then 
X^  -  0.  If  the  flow  is  completely  laminar  on  the  forebody 

surface,  it  may  be  laminar  or  turbulent  in  the  wake.  A  wake 
transition  criterion,  discussed  in  IIB4a(3) ,  roust  be  used  to 
determine  the  location  of  the  transition  point,  X^. 

(3.)  Wake  Transition  Criterion 

The  location  of  the  wake  transition  point  X^  from 

laminar  to  turbulent  wake  flow  can  be  determined  primarily 
by  using  available  experimental  data.  In  Ref  1,  an  approximate 
curve  fit  correlating  experimental  results  of  Pallone,  Erdos, 
and  Eckerman  (Ref  5)  has  been  used  as  the  wake  transition 
criterion.  The  wake  transition  Reynolds  number  Re^^  is 

based  on  viscous  wake  properties  at  the  edge  of  the  wake  with 
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the  relative  velocity  between  the  centerline  velocity  and 
velocity  at  wake  edge  as 


fz  (U,  -  )Jqtr 

yU, 


(47) 


and  effective  Mach  number  is  defined  as 


A1, 


Of,  -  Uct 

a, 


(48) 


At  the  transition  pointy  transition  Reynold's  number  and 
effective  Mach  Number  are  related  (according  to  a  curve 
fit  in  Ref  1)  by 


R 


e 


MOO  M,  2*~ 


(49) 


If  the  flow  on  the  forebody  surface  is  completely  laminar, 
the  laminar  viscous  wake  analysis  of  IIB4a(l)  is  continued 
until  the  wake  transition  criterion  (49)  is  satisfied.  The 
point  XTR  and  the  value  5^aminar  at  this  point  are  used  in 

continuing  the  turbulent  wake  analysis  of  IIB4a(2), 
b.  Inviscid  Wake 


The  inviscid  wake  analysis  essentially  determines 
the  flow  properties  beyond  the  viscous  wake  edge  at  any 
specified  location  behind  the  base  of  the  forebody.  The 
following  assumptions  have  been  made  in  the  far  wake  region. 


1. 


2. 


The  secondary  shocks,  such  as  the  wake 
recompression  shocks,  forebody  flare  shocks, 
if  any.  are  considered  to  be  very  weak  so 
that  they  produce  negligible  change  in  entropy. 


There  is  no  radial  pressure  gradient 
in  the  wake. 


3.  In  the  far  wake  region  static  pressure  in  the 
wake  approaches  the  freestream  static  pressure. 

The  inviscid  wake  depends  mainly  on  the  shape 
and  strength  of  forebody  shock.  The  total  pressure,  total 
temperature  and  entropy  at  any  radius  r  from  the  axis  of 
wake  are  equal  to  the  corresponding  values  at  radius  rfi 

immediately  behind  the  forebody  shock.  Since  the  forebody 


shock  is  completely  specified,  flow  properties  immediately 
behind  the  shock  are  given  by  shock  jump  conditions.  The 
relation  between  r  and  r  is  given  by  the  conservation  of 
mass,  s 


rs*  -  zf  1^-rdr 

o  ^  00 

where  the  integration  in  Eqn  50  is  started  from 
displacement  thickness,  to  take  into  account  the 
of  streamlines.  The  wake  displacement  thickness 

by 


(50) 


r* 

b  ,  the  wake 
displacement 
is  defined 


(51) 


To  facilitate  the  computation  (X,  r)  coordinates  are  trans¬ 
formed  into  (X,  n)  where 


f  Ur  dr 


^oo  ^co  ri  d  n 


(52) 


So  that  Eqn  50  becomes 


2 


fl  U> 

u* 


(53) 


In  order  to  transform  back  into  physical  coordinates,  r  and  n 
are  related  by 


The  computational  procedure  is  as  follows: 


(54) 


1.  From  the  results  of  the  viscous  wake  analysis, 
the  wake  displacement  thickness  &  can  be  calculated  from 
Eqn  51. 

2.  For  any  transformed  radius  n,  Eqn  53  is  used 
to  find  r_. 

a 

3.  Total  pressure,  total  enthalpy  and  entropy  at 

the  transformed  radius  n  are  the  same  as  that  at  radius  r„ 

s 

immediately  behind  the  forebody  shock. 
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4.  Static  pressure  is  assumed  to  be  equal  to 
freestream  static  pressure.  From  total  pressure  and  static 
pressure  the  Mach  number  at  transformed  radius  n  can  be 
found  using  isentropic  relations.  All  other  flow  properties 
at  this  location  can  then  be  found. 

5.  Radial  distance  r  from  the  axis  of  symmetry 
corresponding  to  transformed  radius  n  can  be  found  using 
Eqn  54. 


6.  Since  the  viscous  wake  analysis  depends  on 
the  inviscid  solution  at  the  edge  of  viscous  wake,  an  itera¬ 
tion  procedure  is  adopted  to  determine  the  complete  wake  pro¬ 
file. 


C.  Assumptions  and  Validity 

1.  Forebody  Pressure  Distribution 

The  method  of  predicting  the  forebody  pressure 
distribution  has  been  discussed  in  Appendix  II  of  Ref  1. 

The  surface  pressure  distribution  is  needed  as  an  input  to 
the  boundary  layer  calculations.  The  Newtonian  method  is 
used  if  the  shock  wave  is  detached  from  the  body  and  the 
tangent-cone  method  is  used  if  the  shock  wave  is  attached. 

These  methods  of  predicting  surface  pressure  distribution 
can  be  fairly  accurate  in  many  cases,  but  inaccurate  in 
others.  If  the  error  in  pressure  drag  coefficient  is  chosen 
as  an  indicator  of  the  accuracy  of  these  techniques,  Fig  17.18b 
of  Ref  6  can  be  used  to  estimate  the  range  of  applicability 
of  the  techniques.  From  this  figure  one  finds  that  an  error 
of  about  ten  per  cent  or  less  occurs  if  M^TT  el.  M® 
represents  the  freestream  Mach  number  and  Z  represents  the 
thickness  ratio  of  the  body  under  consideration.  For  the 
Arapaho  C  vehicle,  (Configuration  2.  Fig  2,  Ref  7)  a  thickness 
ratio  of  Z  -  1/8  is  appropriate.  Thus,  we  infer  that  the 
surface  pressure  distribution  technique  of  Ref  1  (for  an 
Arapaho  C  vehicle)  is  valid  only  for  freestream  Mach  numbers 
in  the  range  M®  £  8.  For  the  modified  Arapaho  C  with  a 
flared  afterbody,  tsl/4  (Configuration  3,  Fig  2,  Ref  7), 
the  methods  are  valid  for  M®  *  4.  It  should  be  noted  that 
we  have  used  an  Integrated  pressure  distribution  (pressure 
drag  coefficient)  to  determine  these  ranges  of  validity.  It 
should  be  remembered,  as  pointed  out  in  Ref  6,  page  676,  that 
a  good  agreement  in  pressure  drag  coefficient  does  not  imply 
good  agreement  in  surface  pressure  distribution. 

It  is  difficult  to  estimate  the  errors  of  the 
surface  pressure  prediction  method  in  the  prediction  of  wake 
properties.  However,  in  many  cases  it  would  probably  be 
preferable  to  utilize  the  method  of  characteristics  to 
predict  the  forebody  pressure  distribution. 
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2.  Forebody  Boundary  Layer  Calculations 

The  treatment  of  the  boundary  layer  growth  along 
the  surface  of  the  forebody  is  based  on  empirical  expressions 
and  relationships  that  appear  to  provide  a  good  description 
of  the  experimentally  observed  data.  From  the  standpoint  of 
practical  application,  it  would  appear  that  the  boundary 
layer  techniques  utilized  in  Ref  1  are  adequate  for  engineer¬ 
ing  applications. 

3.  Viscous  Wake  Analysis 

The  assumptions  made  to  arrive  at  a  solution  for 
the  viscous  wake  in  Ref  1  (Fig  97)  are  presented  below: 

1.  The  fluid  is  thermally  and  calorically  perfect, 
and  there  are  no  chemical  reactions  or  other  real 
gas  phenomena. 

2.  Prandtl  number  is  unity  (Pr  =  1). 

3.  Prandtl' s  concept  of  viscous  flow  phenomena 
is  valid  for  the  high-speed  compressible  wake  in 
such  a  way  that  gradients  in  the  streamwise  or 
axial  direction  are  much  smaller  than  those  nor¬ 
mal  to  the  wake  axis,  and  boundary  layer  type 
equations  can  be  used. 

4.  The  details  of  the  base  flow  and  free  shear 
can  be  largely  ignored.  Thus,  the  region  of 
validity  for  the  present  analysis  must  be  considered 
to  extend  from  somewhere  in  the  vicinity  of  the 
wake  neck,  or  further  downstream,  down  into  the 

far  wake. 

5.  Effects  due  to  the  existence  of  an  external 
pressure  gradient  are  negligible. 

These  assumptions  from  the  basis  of  classical 
boundary  layer  theory  have  been  the  basis  for  many  successful 
analytical  studies.  However,  assumptions  3  and  5  become 
increasingly  less  valid  as  tne  base  of  the  forebody  is 
approached. 

A  further  assumption  is  implied  by  Eqn  59  of  Ref  1, 
namely,  that  the  velocity  profiles  in  the  wake  are  similar 
at  any  axial  location.  Tnis  assumption  also  becomes  increas¬ 
ingly  inaccurate  as  the  wake  neck  region  is  approached. 

It  is  also  assumed  in  Ref  1,  Eqn  54,  that  the 
momentum  defect  in  the  wake  is  equal  to  the  momentum  defect 
in  the  boundary  layer  at  the  base  of  the  forebody.  Data 
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presented  in  Table  I  of  Ref  1  shows  that  this  assumption 
is  not  completely  valid.  However,  it  may  be  possible  to 
introduce  an  empirical  factor  that  relates  the  boundary  layer 
momentum  defect  at  the  base  to  the  far  wake  momentum  defect. 

The  choice  of  an  appropriate  turbulent  viscosity 
model  is  discussed  in  Ref  1  and  four  models  are  considered. 

The  third  model, *  k a  S,  («,- uc/)with  k  =  0.005  provided 
the  best  correlation  with  experimental  results.  It  is 
impossible  to  say,  with  only  this  one  body  having  been  con¬ 
sidered,  that  the  same  viscosity  model  will  be  appropriate 
for  all  bodies. 

It  was  found  that  the  viscous  wake  analysis  pro¬ 
vided  agreement  with  experiment  to  well  within  10  per  cent 
in  Ref  1.  Considering  the  number  of  assumptions  employed 
and  the  complex  phenomena  being  studied,  this  is  considered 
to  be  excellent  agreement. 

4.  Inviscid  Wake  Analysis 

The  method  used  to  predict  the  inviscid  wake 
characteristics  in  Ref  1  must  be  considered  only  as  a  first 
order  method.  The  assumption  of  isentropic  flow  along 
streamlines  behind  the  bow  shock  is  valid  only  for  certain 
configurations  where  no  trailing  shock  waves  have  been 
generated.  The  method  of  characteristics  could  be  employed 
for  a  second  order  prediction  of  the  inviscid  wake  region, 
but  this  requires  considerably  more  effort  and  complexity. 

5 .  Summary 

The  forebody  and  wake  analysis  technique,  despite 
some  of  its  shortcomings,  appears  to  be  the  best  available 
method  for  wake  predictions.  Several  improvements  and 
refinements  could  be  made,  but  the  increase  in  accuracy  may 
not  justify  a  more  complext  technique. 

D.  Experimental  Comparisons 

Using  the  methods  just  described,  sample  calculations 
were  completed  for  a  cone-cylinder  body  in  a  supersonic  flow 
with  Mot  »  2.98  and  a  Re/ft  of  8.86  x  105.  The  results  of  these 
calculations  were  checked  with  wind  tunnel  experiments  con¬ 
ducted  at  the  Rosemount  Aeronautical  Laboratories,  University 
of  Minnesota,  at  the  Supersonic  Gasdynamics  Facility  of  the 
Flight  Dynamics  Laboratory  at  Wrighc-Patterson  Air  Force  Base, 
and  at  the  Arnold  Engineering  Development  Center's  von  Karman 
Gas  Dynamics  Facility.  Figure  2  shows  the  geometry  of  the 
forebody  and  the  two  sizes  that  were  tested;  calculations  were 
made  for  the  Db  “  2  in. forebody.  In  this  section 
the  results  of  the  calculations  and  experiments  are  compared. 


20 


1.  Forebody  Shock  Shape 

Although  the  forebody  shock  shape  is  not  a  result, 
it  is  a  required  input  to  the  method  of  calculation.  Thus 
the  foreboay  was  sting  mounted  in  a  Rosemount  Aeronautical 
Laboratory  12  in  x  12  in  blowdown  tunnel  and  Schlieren 
photographs  obtained  (Fig  3).  As  shown,  at  Mo©  =  2,98  the 
shock  is  attached  and  nearly  conical.  From  these  photographs 
the  shock  geometry  was  determined  and  used  as  an  input  to  the 
calculations. 

2.  Forebody  Surface  Pressures 

The  surface  pressures  of  the  forebody  are  a  result 
of  the  calculation.  These  pressures  were  measured  in  tests 
conducted  at  the  AFFDL  Supersonic  Gasdynamics  Facility 
continuous  flow  tunnel  at  Moo  **  2.98.  Figure  4  shows  the 
pressure  distribution  model  installed  in  the  SGF  tunnel, 
and  Fig  3  shows  the  test  results  and  the  results  of  the 
calculation. 

As  shown  in  Fig  5  the  results  are  quite  good.  As 
would  be  expected,  the  measured  surface  pressure  on  the  conical 
portion  of  the  body  agrees  almost  perfectly  with  the  cal¬ 
culated  pressure.  Downstream  of  the  conical  portion  the 
measured  pressures  are  at  first  lower,  then  approach  the 
calculated  values.  This  is  no  doubt  due  to  the  inaccurate 
calculation  of  the  flow  expansion  occurring  at  the  inter¬ 
section  of  the  conical  and  cylindrical  portions. 

3.  Base  Momentum  Defect 

The  momentum  defect  at  the  base  of  the  body  was 
also  measured  with  small  boundary  layer  probes  in  the  tests 
at  the  AFFDL  SGF  tunnel  (Fig  6) .  Figure  7  shows  the  velocity 
profile  results  from  these  measurements  and  Table  1  compares 
the  results  of  these  measurements  with  the  calculation. 

There  is  a  23%  difference  in  these  values,  which  is 
probably  consistent  with  the  accuracy  of  the  assumptions 
made  in  the  boundary  layer  analysis. 

4.  Wake  Characteristics 

Figures  8  through  13  show  comparisons  of  measured 
and  calculated  wake  characteristics.  The  measured  values 
presented  here  have  been  obtained  from  Ref  8.  These  tests 
were  conducted  at  Moo  11  2.98  in  an  Arnold  Engineering 
Development  Center  continuous  flow  tunnel.  One  sees  from 
the  figures  that  there  is  a  considerable  difference  in  the 
results,  due  primarily  to  the  fact  that  the  measured  wake 
is  considerably  thicker  than  the  calculated  wake.  Other 
than  this,  the  general  trends  are  similar  and  the  wake  center- 
line  and  edge  values  are  in  reasonable  agreement.  This 
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Forebody  Pressure  Distribution  Model ,  Sting  Mounted, 


Fig  5  Comparison  of  Measured  and 

Calculated  Forebody  Surface 
Pressures  Mxt2.98  R^a 3362x10s 


Forebody  Base  Momentum  Defect  Measurement 
Boundary  Layer  Probes,  M  =2.98 


0  0.2  0/4  Q6  0.8  1.0 


v/u„ 

Fig  7  Velocity  Profiles  in  the  Boundary 
Layer  of  the  Forebody  at 
its  Base  Mm=298  Rf  ~-  8874xl04 


TABLE  I 

MOMENTUM  DEFECT  AT  THE 
FOREBODY  BASE  AND  IN  THE 
WAKE  ( M,*,  =  2.98,  REYNOLDS  NO/FT 
=  8.855  x  105) 


Location 

Momentum  Defect  ( 11') 

Experimental  Theorv 

■ 

0.081 

0.062 

1/D&  =  ^ 

0.314 

0  .Dt>2 

1/db  ■  7 

0.336 

0.0(2 

i/db  "  9 


0.311 


0.062 
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Fig  8  Wake  Total 
M»  =  2.98 , 


Pressure  Distribution 


Fig  9  Wake  Mach  Number  Distribution 
M»  =  2.98 ,  V[v,  =  5.0 


30 


31 


Fig  11  Wake  Mach  Number  Distribution 
Moo  =2.98,  Vdo'7*0 
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Fig  12  \Afeke  Total  Pressure  Distribution 
Moo  =  2.98,  i/D  =9.0 


—  Experimental 
Result  (Ref  8) 


Fjg  13  Wake  Mach  Number  Distribution 
Moo  =  2.98,  ?/n  =9.0 


m '■ 
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discrepancy  in  wake  width  could  be  lessened  by  adjusting 
some  of  the  coefficients  used  in  the  wake  analysis,  but  time 
did  not  allow  this  effort. 
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III.  SECONDARY  BODY  ANALYSIS 


In  view  of  the  over-all  objectives  of  this  study, 
once  the  analysis  of  the  primary  body  had  been  completed, 
the  next  effort  was  directed  at  the  secondary  body,  the 
decelerator.  Analytical  techniques  have  been  developed  for 
predicting  the  flow  field  properties  about  certain  classes 
of  secondary  bodies  using  the  method  of  characteristics. 

This  section  explains  the  theory  of  the  method  used  and 
compares  its  results  with  measurements;  details  of  the 
actual  calculations  are  described  in  Vol  II  of  this  report. 

A.  Extent  of  the  Secondary  Body  Analysis, 

Required  Inputs,  and  Results 

1.  Extent  of  the  Analysis 

The  analysis  used  for  the  secondary  body  extends 
axially  from  the  upstream  tip  of  the  secondary  body  down¬ 
stream  to  the  base  of  the  body,  The  wake  of  the  secondary 
body  is  not  included  in  the  analysis.  Normal  to  the  axis  of 
symmetry,  the  analysis  extends  from  the  secondary  body 
surface  to  its  bow  shock  wave;  flow  field  characteristics 
outside  of  the  bow  shock  can  be  obtained  from  the  wake  analysis 
of  the  forebody. 

2.  Required  Inputs 

The  information  required  for  the  secondary  body 
analysis  includes  the  upstream  flow  properties,  the  pressure 
and  temperature  distributions,  and  the  geometry  of  the  body. 

The  computer  program  discussed  in  Vol  II  is  such  that  the 
upstream  flow  conditions  can  either  be  calculated  using  the 
primary  body  analysis  or  an  input  to  the  program.  In  order 
to  obtain  a  pressure  drag  coefficient  the  base  pressure  of 
the  body  must  be  obtained, 

3.  Results  of  the  Analysis 

The  method  calculates  all  the  flow  properties  from 
the  body  surface  to  the  shock  wave.  Thus  the  shock  shape, 
surface  pressures,  and  pressure  drag  are  obtained. 

B.  Explanation  of  the  Method  Used  for  the  Secondary 

Body  Analysis 

1.  Theory 

In  order  to  predict  the  flow  field  surrounding  an 
aerodynamic  decelerator  in  the  wake  of  a  forebody,  several 


assumptions  and  restrictions  have  been  made. 

1.  The  flow  is  inviscid  and  steady 

2.  The  gas  is  thermally  and  calorically  perfect. 

3.  The  flow  is  supersonic  everywhere  in  the 
region  of  interest. 

4.  The  wake  flow  field  is  known,  either  from 
theory  or  experiment. 

The  first  three  assumptions  permit  one  to  develop 
a  method  of  characteristics  for  rotational,  non-homoenergic 
flow.  The  fourth  assumption  establishes  an  effective  "up¬ 
stream"’  boundary  condition  for  the  problem. 

The  basic  equations  under  the  above  assumptions 
are  as  follows: 

Continuity  equation, 

V-  tV  =  O  (55) 


Momentum  equation, 

T(v-v)v 


-VP 


First  Law  of  Thermodynamics, 

vh  *  TV  £  + 


Using  the  vector  identity 


(V-  V)V  *  -V  x(V  X  V  )  +  7  ~ 

4* 


Equation  56  can  be  written  as 

1  f 

where  co  represents  the  vorticity  and  is  defined  by 


OO  -  V  X  V 
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Solving  Eqn  57  for  and  substituting  into  Eqn  59  yields 
a  convenient  form  of  the  momentum  equation 

Vxa?  =  -TVS  +  v(h+  y)  (61) 

If  one  takes  the  dot  product  of  the  velocity  \/  with  Eqn  61, 
there  results 


-T  V -V  S  *  V'V(h+  j_V2)  --  O  (62) 

If  one  assumes  that  the  total  energy  is  constant  along  a 
streamline,  then  Eqn  62  shows  that  the  entropy  is  also 
constant  along  a  streamline. 

Taking  the  dot  product  of  Eqn  61  with  a  unit  nor¬ 
mal  vector  to  the  streamline,  -A  ,  one  obtains 


n-(vx<2)  *  -T 


65  *  ih* 


where  h^  ==  h  +  ^  “  total  enthalpy.  (( 

The  vorticity,  in  cylindrical  coordinates,  can 
be  written  as 


It  can  be  shown  that,  for  axially  symmetric  flow,  the  unit 
normal  vector  is  given  by 


■f 


V 


(66) 


One  can  then  evaluate  the  left  sideaf  Eqn  63  to  obtain 

\  h*.  -  lit  ]  W  r  -  T  ^5.  +  4-b®. 

lar  v  d*  dr\  (67) 
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Using  conventional  relationships  between  local  static  and 
total  properties  (Ref  2),  it  can  be  shown  that  the  density 
and  speed  of  sound  are  given  by 


J _ 

2M 


(68) 


(69) 


Substituting  Eqn  68  into  the  continuity  equation  55  gives 

i 


-p. 


(2.  h e )  y-i 


(2k  -vl) 


which  may  be  expanded  to 


(70) 


(2  k  -V')  U 

V-[(2k-Vl) 


(71) 


Since  total  pressure  (pQ)  and  enthalpy  (hQ)  have  been  shown 
to  be  constant  along  a  streamline,  one  can  write 


V 

and  Eqn  71  becomes 

V 


2 

(2  h,  )  * 


=  O 


=  o 


(72) 


(73) 
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In  cylindrical  coordinates,  this  can  be  expressed  as 


[y  (zKr^) 


(74) 


This  relation  implies  the  existence  of  a  stream 
function  that  satisfies  the  relations 


A  characteristic  equation,  involving  the  stream 
function,  can  be  developed  from  Eqn  67  in  the  following  fashion. 
As  a  convenient  shorthand  notation,  we  will  indicate  partial 
differentiation  of  a  variable  with  a  subscript.  Thus,  for 
example, 


l Y  _ 


t± 

dxdr 


(77) 

(78) 


If  we  let 


[*k»  - 


-c 


(79) 


Eqns  75  and  76  can  be  written 

U  =  j  fr 


(80) 


*  -f  % 


(81) 
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Using  these  relations,  the  bracketed  portion  on  the  left 
side  of  Eqn  67  can  be  written 


ur  -  ^  =  (%  *  frr)j  +(jrfr*  f,  %)  (82) 

By  squaring  Eqns  80  and  81  and  adding  the  results,  one 
obtains 

yC  ■>  =  ji  <83) 

Taking  the  gradient  of  Eqn  83  gives 

y* 

-  v( Ji)  (84) 


Taking  the  dot  product  of  n,  Eqn  66  with  Eqn  84  results  in 


t. 


- 

V  h 


f* 


X/* 


U * 

f* 


One  can  rearrange  Eqn  79  to  obtain 


,  -(**<) 
r  f 


(85) 


(86) 


Using  this  relation,  the  right  side  of  Eqn  80  can  be  evaluated. 
After  considerable  algebraic  manipulations,  Eqn  85  is  trans¬ 
formed  to  the  following 


fxf*  *  frft 
•f 


-  +  tuvfyir  -  “*fy, 

-t  ~?*i 
f  V 


rr 


tJ.Y' 
%  T 


<7rT 


(87) 
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Substituting  Eqn  87  into  Eqn  78  provides,  after  algebraic 
manipulation 


»  j^(u2-  o?)  ^  +  ('isl-a2)  frr  +  zu^fxr  + 

+  ~  r  /  'tc.  +  V~  d  In  o']  f 

r  f  dn\  Vl-al 


Substitution  of  Eqn  88  into  Eqn  67  provides  the  final 
equation 


Z 


U.V~ 

a1 


5?  ft*' i 


t-$l 


d  % 

dm 


(89) 


The  equation  for  the  stream  function  is  identical  to  Eqn  9.65 
of  Ref  6  with  the  exception  of  the  terms  on  the  right  side 
of  Eqn  89.  Equation  9.65  of  Ref  6  was  derived  for  irrotational 
flow.  Thus,  the  terms  on  the  right  side  of  Eqn  89  describe 
the  effect  of  entropy  and  enthalpy  gradients  on  the  flow  field. 

The  general  theory  of  characteristics  is  presented 
and  discussed  in  Appendix  A  of  Ref  6,  A  characteristic  curve 
is  one  on  which  and  %  are  discontinuous.  If  we  write 
Eqn  39  in  the  form 


A %  +  28  ftr  +  C.  fn.  -  °  O0) 

it  can  be  shown  that  the  slope,  "X  ,  of  the  characteristic 
curves  is  given  by 

\  -  dr  8-  AC  (91) 

"  - A 
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As  shown  in  Ref  6,  page  484,  this  relation  reduces  to 


>=  =  ian  (&*/<)  («) 

where yuc  is  the  Mach  angle  given  by 

sm  ju.  x  -L  (93) 

/  M 

Since  siri/^.  is  defined  only  for  M  2  1,  one  notes  an  important 
feature  of  the  method  of  characteristics,  namely,  it  can  be 
used  only  in  regions  where  the  flow  field  is  supersonic. 

We  now  proceed  to  establish  the  equation  describing 
the  variation  of  the  stream  function  along  the  characteristic 
curves.  In  general,  one  may  write 


d\  ’  o'*  +  -  1% +  f.r  fc]dX 


(94) 


On  a  characteristic  curve,  Eqn  94  can  be  written 


(&L-  ^ 


J char 

which  can  be  rearranged  to  give 


(95) 


*  /(Mar 

In  a  similar  manner,  one  finds 


(96) 


If  one  substitutes  Eqns  96  and  97  into  Eqn  90  and  performs 
some  algebraic  manipulations,  the  following  relations  can 
be  derived*. 


on  ; 

d% 


-dj  +  1  d% 


D  _ 


—  Q 


dr.  =  ;  dv,  A.  djr  _  A  =  O 

c/%  o'*.  ^ 


Using  Eqn  75  one  can  write 


2\ 


yr  =  r~(2b0-V  )  V  cos  G 


(100) 


Differentiation  yields 


5T-/^7: 


-  r  4&  +  r  cos  &  4v.  + 

dx 


+  ir  ~^a  ihtek-v')  +  COS  &] 

t-l  2h0-v7dx[  /  j 


(101) 


Similarly  one  obtains 


--(2A.-0  *"f- 


Y  V  Cose  4&  ~  r  str)  &  zr 
dx  dx 


r->  tK.„*  d^zh° v> 


A  Vs  »«e 


(102) 
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If  one  substitutes  Eqns  101  and  102  into  Eqns  98  and  99,  it 
can  be  shown,  after  a  considerable  amount  of  algebra,  that 
the  equations  along  the  characteristics  can  be  expressed  as 


* 

-tan  (  e  dyu.) 


(103) 
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Introducing  D  and  A  from  Eqn  89  (by  comparison  with  Eqn  90), 
one  can  obtain  the  final  characteristic  equations. 
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The  upper  sign  in  Eqns  105  and  106  refer  to  the  so-called 
left-running  characteristic  while  the  lower  sign  refers  to 
the  right-running  characteristic. 


If  one  multiplies  Eqn  106  by  dx  and  notes  that 


Equation  106  becomes: 
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Eqn  108  is  integrated  along  a  characteristic 
to  Point  2,  and  if  the  average  value  of  the 
for  each  term  are  used,  Eqn  108  becomes 
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All  quantities  with  a  "bar"  over  them  and  subscript  12 
represent  average  values  between  Points  1  and  2  on  a 
characteristic  curve.  Once  again  it  should  be  remembered 
that  the  upper  signs  are  used  for  the  left  running 
characteristics  while  the  lower  signs  are  for  the  right- 
running  characteristics. 

2.  Computational  Procedure 

A  detailed  discussion  of  the  computation  in  the 
shock  layer  flow  field,  such  as  the  location  of  shock  points. 
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net  points,  body  points,  and  the  properties  along  an  initial 
right-running  characteristic  is  presented  in  Vol  II  of  this 
report.  A  brief  description  of  the  flow  field  computation 
has  been  included  in  this  section  to  provide  a  better  under¬ 
standing  of  the  characteristic  and  compatability  equations. 

As  described  earlier,  the  non-uniform  supersonic  wake  flow 
field,  which  is  known  either  from  an  analytical  approach 
(Section  II)  or  from  experiments,  provides  effective  upstream 
boundary  conditions;  the  Rankineflugonoit  relations,  when 
applied  locally  across  the  shock,  determine  the  flow  proper¬ 
ties  immediately  behind  the  shock.  With  an  attached  curved 
shock  at  the  apex  of  the  secondary  body,  the  shock  layer, 
the  region  between  the  shock  and  body  surface,  will  have 
rotational,  non-homoenergic  characteristics.  The  analysis 
using  the  method  of  characteristics  is  valid  in  the  shock 
layer  provided  the  flow  field  is  completely  supersonic. 

In  order  to  proceed  with  the  computation  of  the 
flow  field  by  the  method  of  characteristics,  it  is  necessary 
to  know  the  flow  properties  along  an  initial  right  running 
characteristic.  In  the  case  of  an  attached  shock,  it  is 
usually  assumed  that  a  small  conical  flow  region  exists  at 
the  apex  of  the  body.  The  shock  angle,  related  to  the  semi¬ 
apex  angle  of  the  body  and  the  centerline  upstream  Mach 
number,  and  the  flow  properties  in  the  conical  flow  region 
are  determined  from  Taylor-Maccoll  conical  flow  relationships. 
Since  the  local  characteristic  direction  is  inclined  at  the 
local  Mach  angle  from  the  local  flow  direction,  the  initial 
right  running  characteristic  can  be  constructed  from  the 
knowledge  of  the  local  Mach  number  and  local  flow  direction 
in  the  conical  flow  region.  The  computation  along  any  right 
running  characteristic  is  carried  out  from  a  shock  point  to 
a  body  point.  Knowing  the  location  and  flow  properties  at  a 
shock  point,  net  points,  and  a  body  point  on  one  right  running 
characteristic,  the  location  and  flow  properties  at  field 
points  on  the  next  right  running  characteristics  are  deter¬ 
mined  using  shock  jump  conditions,  characteristic  compatability 
conditions,  and  body  boundary  conditions. 

A  new  shock  point  is  fixed  by  using  the  shock  point 
and  the  net  point  on  the  previous  right  running  character¬ 
istic  and  Rankine-Hugonoit  shock  jump  conditions.  An  itera¬ 
tion  procedure  has  to  be  adapted  to  fix  the  new  shock  point. 

A  new  net  point  is  located  using  the  previously 
located  shock  point  and  net  point,  two  net  points,  or  net 
point  and  body  point,  depending  on  the  region  of  desired  new 
net  point.  In  addition,  the  streamline  passing  through 
the  net  point  in  search  must  be  located  on  the  previous 
right  running  characteristic  to  determine  the  total  tempera¬ 
ture  and  entropy  at  the  new  net  point.  The  velocity  and  flow 
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direction  at  the  new  net  point  are  determined  from  two  char¬ 
acteristic  compatability  relations  (Eqn  109) . 

Since  the  body  surface  is  a  streamline,  the  total 
temperature  and  entropy  remain  constant  along  the  body  surface. 
From  body  boundary  conditions  the  flow  direction  at  any  point 
on  the  body  surface  is  along  the  tangent  direction  to  the 
surface  at  that  point.  The  new  body  point  is  located  using 
the  previously  determined  body  point  and  net  point.  The 
compatability  relation  along  the  right  running  characteristic 
provides  the  magnitude  of  the  surface  velocity.  The  tempera¬ 
ture  and  static  pressure  are  determined  from  the  isentropic 
relations 
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Then  the  local  velocity,  flow  direction,  total  temperature 
and  entropy  are  sufficient  to  determine  all  other  flow  pro¬ 
perties  at  any  point. 


C.  Limits  of  and  Assumptions  in  the  Secondary  Body 

Analysis 

1.  Limits  of  the  Secondary  Body  Analysis 

The  primary  and  obvious  limitation  is  that,  since 
a  characteristic  metnod  is  used,  it  is  limited  to  supersonic 
velocities.  The  most  important  implications  of  this  fact 
are  that  the  secondary  body  must  have  an  attached  bow  shock 
and  the  analysis  will  only  extend  along  the  body  to  the 
sonic  line.  Although  this  limitation  is,  in  a  certain  sense, 
severe,  there  is  yet  a  large  number  of  bodies  to  which  the 
method  will  apply. 

A  second  limitation  of  the  method  is  that  the 
wake  of  the  secondary  body  is  not  analyzed.  The  result  of 
this  is  that  base  pressures  must  be  obtained  from  another 
source  to  obtain  a  complete  pressure  distribution  for  a  drag 
coefficient  calculation. 
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2.  Assumptions  in  the  Method 

Two  main  assumptions  are  made  in  using  the  method: 
first,  it  is  assumed  that  a  boundary  layer  analysis  is  not 
needed  to  obtain  the  desired  results,  and  second,  that  a 
reasonably  sized  decelerator  tow  line  will  not  invalidate  the 
calculation  with  an  attached  shock. 

The  omission  of  a  boundary  layer  analysis  should 
not  lead  to  any  significant  error  in  the  results  obtained. 

For  bodies  and  regimes  to  which  the  method  applies,  boundary 
layer  thicknesses  are  quite  small  compared  to  body  dimensions, 
and  the  friction  drag  is  a  small  percentage  of  pressure  drag. 

The  primary  need  for  a  boundary  layer  analysis  would  be  for 
wake  calculations. 

At  the  very  early  stages  of  this  investigation 
wind  tunnel  studies  were  conducted  to  examine  the  effects 
of  simulated  axial  tow  lines  on  the  secondary  body  bow  shock 
shape  and  pressure  distribution.  Figure  14  compares  the 
shapes  of  bow  shock  waves  on  a  4-in  diameter  30°  half-angle 
cone  in  the  wake  of  the  2-in  diameter  cone-cylinder  forebody 
with  and  without  two  different  simulated  tow  lines.  As  shown 
in  Fig  14  the  bow  shocks  with  the  simulated  tow  lines  are 
obviously  not  attached,  but  the  region  of  difference  is  very 
small,  and  the  shock  rapidly  approaches  that  of  the  no-tow- 
line  configuration.  Pressure  distribution  studies  were  also 
conducted  for  this  same  cone  with  and  without  a  tow  line,  and 
the  results  are  shown  in  Fig  15.  It  can  be  seen  that  the 
simulated  tow  lines  have  nearly  no  effect  on  the  pressure 
distribution. 

It  should  be  noted  that  one  effect  of  the  simulated 
tow  line  was  not  directly  investigated.  This  effect  is  that, 
for  a  thickness  equal  to  that  of  the  tow  line,  there  exists 
a  wake  velocity  boundary  condition  of  V  ■  0.  This  was  not 
considered  in  either  the  wake  analysis  or  the  secondary  body 
analysis.  This  condition  will  exist  in  any  towed  decelerator 
application,  but  the  effects  of  neglecting  it  cannot  be  iso¬ 
lated  in  the  results  of  this  study. 

D.  Experimental  Comparisons 

In  order  to  check  the  results  of  a  sample  calculation 
using  the  analysis,  experiments  were  conducted  to  measure 
secondary  body  bow  shock  shape  and  surface  pressure  distribu¬ 
tion.  These  tests  were  conducted  at  the  University  of  Minnesota's 
Rosemount  Aeronautical  Laboratories  in  a  12  in  x  12  in  blow¬ 
down  wind  tunnel  at  M*  -  3.0.  Figure  16  is  a  schematic 
of  the  installation  and  Fig  17  shows  the  30°  half-angle  cone 
used  as  the  secondary  body.  The  primary  body  had  the  same 
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Fig  14  Comparison  of  Bow  Shock  Wave 
Shapes  with  and  without  two 
Simulated  Towline  Diameters  at 
fi /  Dg  =  5 1  M*  =2.98 
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Fig  15  Pressure  Distribution  on  a  30° 

Half-Angle  Cone  J/D*5,  NV2.98 
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geometry  as  the  1.2-in  diameter  cone-cylinder  body  examined 
in  Section  II,  and  the  tests  were  conducted  at  the  same  condi¬ 
tions.  Measured  shock  shapes  and  surface  pressure  distribu¬ 
tions  were  non-dimensionalized  with  the  2. 4- in  diameter  cone 
and  converted  to  shapes  and  distributions  for  the  4 .0-in 
diameter  cone.  This  was  done  to  allow  direct  comparison  with 
the  theoretical  results  for  the  4-in  diameter  cone  and  2-in 
base  diameter  cone-cylinder  forebody. 

Generally,  in  the  comparisons  that  follow,  two 
different  calculations  of  the  various  parameters  were  made. 

The  first  calculation  was  made  using  wake  data  obtained  from 
measurements  (Ref  8) ,  and  the  second,  used  wake  data  resulting 
from  the  forebody  wake  analysis.  This  was  done  in  order  to 
provide  a  check  of  the  secondary  body  analysis  in  which  all 
the  inputs  to  the  analysis  are  as  accurate  as  possible.  It 
is  shown  in  the  following  sections  that  better  agreement  is 
obtained  using  measured  rather  than  calculated  wake  data. 

1.  Secondary  Body  Bow  Shock  Shape 

Shock  shape  measurements  were  made  at  ^/Dg's  of 
5,  7,  and  9  (Figs  18,  19,  20).  At  an  //Dg  of  5  experimental 
data  was  not  available,  so  the  comparison  only  compares  measured 
and  calculated  points  using  the  primary  body  analysis.  The 
agreement  at  J/Dg’s  of  7  and  9  between  measured  values  and 
values  using  measured  wake  data  is  excellent,  very  nearly 
exact.  The  agreement  with  the  purely  calculated  shock  shape 
is  not  as  close,  but  must  be  called  good. 

2.  Surface  Pressure  Distribution 

.  Pressure  distributions  were  also  measured  at 

x/Dg's  of  5.  7,  and  9,  as  shown  in  Figs  21,  22,  and  23. 

Again  the  agreement  is  best  with  the  calculations  using  the 
experimental  wake  data,  being  within  101  or  better.  The 
effects  of  the  calculated  wake  being  too  narrow  are  seen  in 
the  Cp  results  based  on  this,  in  that  the  Cp’s  rise  too 
rapidly  from  the  lower  values  in  the  front  of  the  cone. 

It  appears  that,  if  the  width  of  the  predicted  wake  was 
adjusted,  the  agreement  could  be  brought  into  the  10%  region. 

3.  Pressure  Drag  Coefficients 

Calculated  and  measured  drag  coefficients  are 

shown  in  Fig  24.  In  general  the  calculated  drag  coefficients 
are  11%  higher  than  measured,  and  the  difference  between 
values  using  measured  or  calculated  wake  conditions  is  quite 
small. 


A  measured  base  pressure  coefficient  of  Cn  »  -0.074 
was  used  to  obtain  the  pressure  drag  coefficients.  ' 
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Fig  20  Shock  Shapes  fora  30°Ha!f-Angle  Gone  in  the  Wake 
of  the  Forebody  at  =9.0,  =  2.98 
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Fig  23  Surface  Pressure  Coefficient  vs  Axial  Distance 

for  a  30°  Half  -  Angle  Cone  Behind  a  26°  Cone- 
Cylinder  1  /  DR  =  9.0  ,  Mffi  =  2.98 
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Fig  24  Pressure  Drag  Coefficient  vs  I  /  DR 
for  30°  Half- Angle  Cone  Behind  a 
26°  Cone-Cylinder .  Moo*  2.98 


IV.  SUMMARY 


Methods  of  analysis  have  been  developed  for  calcu¬ 
lating  flow  field  properties  of  certain  classes  of  two-body 
decelerator  systems  in  supersonic  flow.  In  this  volume  the 
analytical  techniques  used  for  the  forebody  and  secondary 
bodies  have  been  explained,  and  the  results  of  a  sample  cal¬ 
culation  compared  with  experimental  measurements.  In  general, 
the  comparison  was  quite  good;  differences  between  calculated 
and  measured  values  were  on  the  order  of  10%. 

Thus  methods  and  procedures  exist  for  calculating 
the  significant  flow  field  properties  of  the  two  body  system 
with  reasonable  engineering  accuracy.  In  achieving  this 
the  information  about  the  flow  field  required  for  the  cal¬ 
culations  has  been  kept  to  the  smallest  amount  possible  at 
this  time.  However,  two  things  are  needed  in  order  to 
calculate,  namely,  the  forebody  bow  shock  shape  and  the 
secondary  body  base  pressure.  These  two  items  can  easily 
be  determined  experimentally,  or  can  be  approximated  quite 
well  for  initial  calculations. 

In  conclusion,  a  large  step  has  been  made  in 
enabling  calculation  of  the  performance  of  certain  supersonic 
decelerators,.  Certain  improvements  and  refinements  can 
and  should  be  made,  but  the  broad  framework  and  basic  calcu¬ 
lations  have  been  established  making  modifications  possible. 
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